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Abstract 

We investigate the generalized parton distributions (GPDs) with non-zero ^ and A-*- for a rel- 
ativistic spin-1/2 composite system, namely for an electron dressed with a photon, in light-front 
framework by expressing them in terms of overlaps of light-cone wave functions. The wave function 
provides a template for the quark spin-one diquark structure of the valence light cone wave function 
of the proton. We verify the inequalities among the GPDs with different helicities and show the 
qualitative behaviour of the fermion and gauge boson GPDs in the impact parameter space. 
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I. INTRODUCTION 



Generalized parton distributions (GPDs) have attracted a considerable amount of theo- 
retical and experimental attention recently. An interesting physical interpretation of GPDs 
has been obtained in 0, I2I by taking their Fourier transform with respect to the trans- 
verse momentum transfer. When the longitudinal momentum transfer ^ = 0, this gives the 
distribution of partons in the nucleon in the transverse plane. They are called impact param- 
eter dependent parton distributions (ipdpdfs) q{x,b^). In fact they obey certain positivity 
constraints which justify their physical interpretation as probability densities. This inter- 
pretation holds in the infinite momentum frame (even the forward pdfs have a probabilistic 
interpretation only in this frame) and there is no relativistic correction to this identification 
because in light-front formalism, as well as in the infinite momentum frame, the transverse 
boosts act like non-relativistic Galilean boosts. It is to be remembered that the GPDs, 
being off-forward matrix elements of light-front bilocal currents do not have a probabilistic 
interpretation, rather they have interpretation as probability amplitudes. q{x, 6"*") is defined 
in a proton state with a sharp plus momentum p"*" and localized in the transverse plane such 
that the transverse center of momentum = (normally, one should work with a wave 
packet state which is very localized in transverse position space, in order to avoid the state 
to be normalized to a delta function . q{x,b^) gives simultaneous information about 

the longitudinal momentum fraction x and the transverse distance b of the parton from the 
center of the proton and thus gives a new insight to the internal structure of the proton. 
The impact parameter space representation has also been extended to the spin-dependent 
GPDs [l| and chiral odd ones 4 1 . 

GPDs Hg{x,0,t) have been investigated in the impact parameter space in several ap- 
proaches, for example in the transverse lattice formalism for the pion , in a two component 
(spectator) model ^ for the nucleon, in the chiral quark model for the pion ,7] and using a 
power law wave function for the pion P]. The spin- flip GPD Eg has not been addressed in 
these. The connection of Eg in the impact parameter space and the Siver's effect has been 
shown in P| within the framework of the scalar diquark model of the proton. In a previous 
work [1^, we have calculated both H{x,0,t) and E{x,0,t) in the impact parameter space 
for a spin-1/2 composite relativistic system, namely for an electron dressed with a photon in 
QED. The state can be expanded in Fock space in terms of light-cone wave functions. The 
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GPDs are expressed as overlaps of light-cone wave functions The wave functions in this 
case can be obtained from perturbation theory, and thus their correlations are known at a 
certain order in the coupling constant. Their general form provides a template for the effec- 
ive quark spin-one diquark structure of the valence light-cone wave function of the proton 



12j |. Such a model is self consistent and has been used to investigate the helicity structure 
of a composite relativistic system An interesting advantage is that the two-body Fock 
component contains a gauge boson as one of its constituents and so it is possible to inves- 
tigate the gauge boson GPDs Hg and Eg. Studies of the deep inelastic scattering structure 
functions in this approach and for a dressed quark state also yield interesting results [3, 1^ • 

So far we have discussed GPDs in impact parameter space for ^ = 0. However, deeply 
virtual Compton scattering experiments probe GPDs at nonzero ^. In this case, a Fourier 
transform with respect to the transverse momentum transfer A"*- is not enough to diagonalize 
the GPDs and thus giving a density interpretation. As the longitudinal momentum in the 
final state is different from that in the initial state, the resulting matrix element would still 
be off-diagonal. Recently, certain reduced Wigner distributions, when integrated over the 
transverse momenta of the partons are shown to be the Fourier transforms of GPDs and 
they can be interpreted as the 3D density in the rest frame of the proton for the quarks with 
light cone momentum fraction x. In fact, integration over the z coordinate relates them to 
the ipdpdfs with ^ = 0. In P] it has been shown that for nonzero ^, the Fourier transform of 
the GPDs with respect to probes partons at transverse position 6^, with the initial and 
final protons localized around 0^ but shifted from each other by an amount of order ^6-*-. At 
the same time, the longitudinal momentum of the protons are specified. This difference of 
the transverse position of the protons depends on ^ but not on x and thus this information 
should be present in the scattering amplitudes measurable in experiments where the GPDs 
enter through a convolution in x. This aspect makes it interesting to investigate the GPDs 
in the impact parameter space for nonzero ^. Here also, a useful approach is based on 
the overlap representation of GPDs in terms of light-cone wave functions The overlap 
representation can also be formulated directly in the impact parameter space, in terms of 
overlaps of light-cone wave functions ip{x, b^), which are the Fourier transforms of the wave 
functions with definite transverse momenta ipi^j k-^)- 

Here, we calculate the GPDs Hg g{x,^,t) and Eq g{x,^,t) for an effective spin-1/2 system 
of an electron dressed with a photon in QED and we investigate them in the impact param- 



3 



eter space. The plan of the paper is as follows. The definitions of the fermion and gauge 
boson GPDs are given in section II. The fermion and gauge boson GPDs are calculated 
respectively in section III and IV for a dressed electron state. The GPDs are expressed in 
the impact parameter space in section V. The issue of certain inequalities among the GPDs 
in the impact parameter space is addressed in section VI. The summary and discussions are 
given in section VII. 



II. GENERALIZED PARTON DISTRIBUTIONS 



The GPDs are defined in terms of off-forward matrix elements of light-front bilocal cur- 
rents. In the light-front gauge A'^ = we have, 
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where the ellipses indicate higher twist terms. The momenta of the initial (final) state is 
P{P') and helicity A(A'). Ux{P) is the light-front spinor for the proton. The momentum 
transfer is given by A'^ = P'^' - P^, skewedness ^ = - At. The average momentum of the 
initial and final state proton is P^ 



E!l±£l!L_ Wg take the frame where P-*- = 0. Without 



any loss of generality, we take ^ > 0. t is the invariant momentum transfer in the process, 

t = A^. For simplicity we suppress the flavor indices. Following ^ we define 

p/± p± ^± 



D- 



(2.3) 
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Hng and Eq^g are the twist two fermion and gauge boson GPDs. Using the light-cone spinors 
we get 

e 



2M7r^F 



Eq^g{x,i,t). 



(2.4) 



(2.5) 



Note that Eq^g appear both in helicity-fiip and helicity non-flip parts. For ^ = 0, if, 
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correspond to nucleon helicity non-flip and Eq „ correspond to the helicity flip part 



III. FERMION GPDS 



We take the state | P, a) to be a dressed electron consisting of bare states of an electron 
and an electron plus a photon : 

|P,a) =0i6t(p,^)|o) 

MP, ^ I h, (Ti; k2, X2p{ki, cri)a\k2, X2) \ 0). (3.1) 

Here and are bare photon and electron creation operators respectively and 0i and 
{/)2 are the multiparton wave functions. They are the probability amplitudes to find one bare 
electron and one electron plus photon inside the dressed electron state respectively. 

We introduce Jacobi momenta Xi,qi^ such that J^i^i = 1 and J^iQi^ = 0. They are 
defined as 

x, = -^, q^ = k^-x.P^. (3.2) 
Also, we introduce the wave functions, 

A = 01, Mxi,q^) = ^HKM^)\ (3.3) 

which are independent of the total transverse momentum P-*- of the state and are boost 
invariant. The state is normalized as, 

(P', A' I P, A) = 2{2TifP+6x,x'5{P+ - P'^)6\P^ - P'^). (3.4) 

The two particle wave function depends on the helicities of the electron and photon. Using 
the eigenvalue equation for the light-cone Hamiltonian, this can be written as [iS^ . 

/ In Xil — x) 1 e \ \r. ■ ~ \ 



(g^)2 + m2(l - xf ^(l-x) y^2(27r)3 
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xA*^i- (3.5) 
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m is the bare mass of the electron, 5"^ = —a^ and = a"^. ipi actually gives the normaliza- 



tion of the state 



t/)! 1^ = 1 - / dx \ ' log — - ^ + (1 - a;) , (3.6) 



OL r^-' , rl + x"^ r 1 + 



2n 



1 — x m'^{l — x)'^ 1 — x 



within order a. Here e is a small cutoff on x, the longitudinal momentum fraction carried 
*e fe^ion. We Kave take, the cutoff of the t.a„sve.e .o,„e.ta to be Q. This 
gives the large scale of the process. The above expression is derived using Eqs (|3.4p . (j3.ip 
and (jS3D. 

The helicity non-flip part of the matrix element F^X gives information about both 
and Eq^ as can be seen from ()2.4|) . In terms of the wave function this can be written as, 



= I 1^5(1 -x) + J dVV^;^(^, + (1 - x)D^)^,^{^, q^) 



■ + (3-7) 

We restrict ourselves to the DGLAP region 1 > x > ^. As we are considering no antiparticles 
< a; < 1 in our case. It is known that in the ERBL region, — ^ < a; < ^ the GPDs 
are expressed as off-diagonal overlaps of light-cone wave functions involving higher Fock 
components 

The normalization of the state, | ipi |^ given by eq. ()3.(ij) gives another 0{a) contribution. 
The g-*" integral in the above expression is divergent and can be performed using the same 
cutoffs as discussed above. We get, using Eq. ()3.5p . 
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A2 

~ ''^''^^(l-x)2Di + m2(l-a;i)2] 
+ 2m2(l -xi)2(l-X2)'/}, 



\J{1- Xl){l - X2) 

m^{l-X2fl 



{2(1 + 



X1X2) 



'I + X1X2) Tclog[- 



A^ 



m^(l — 0:2)^ 
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(3.J 



where I = J f^, Li = (g^ + (1 - + m2(l - xif and L2 = (q^f + m\l - X2) 



2. 



Xi — X2 
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It is especially interesting to investigate ()3.7p in the forward limit. For simplicity, we 
consider the massless case. We get. 



F+|(x,0,0) = \ iJi f6{l-x)+ j d'q^r2+{x,q^)i^2+{x,q^) 



2 a 1 + x^ 
' 27r 1 -X 



log-,, 



(3.9) 



here is a scale, /i << A. 

The normalization in this case gives, 

1^1^ = 



a 
2^ 



1-^ , 1 + x2 A2 
ax log—. 

/i2 
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(3.10) 



Thus we have, 

F+^(x,0,0) = 5{l-x) + -log-[j^^^^+-5{l-x)\. (3.11) 

Here the plus prescription is defined in the usual way. Now we know that in the forward 
limit, 0,0) = q{x) which is the (unpolarized) quark distribution of a given fiavor in 

the proton. From ()2.4p we then get F^l{x, 0, 0) = q{x), so from ()3.11|) we get the splitting 
function for the leading order evolution of the fermion distribution 

P,,{x) = (3.12) 

Note that the e dependence is no longer there in Eq. 1)3.111) . This result is also obtained in 
1^ by calculating the structure function of a quark dressed with a gluon (here one would 
also get the color factor Cf). In the nonforward case, one obtains the splitting function for 
the LO evolution of the GPDs Finally, it can be shown from (13. lip that 

r dxHg{x, 0, 0) = Fi(0) = 1, (3.13) 

^ 

where Fi{0) is the Dirac form factor at zero momentum transfer. 
Next we calculate the helicity fiip matrix element 

fle^^"^- {P'+ I i^i-^-h^U I P-). (3.14) 

Contribution to ()3.14|) comes from the two-particle sector of the state. The mass cannot 
be neglected here. It can be written as 

- f t/V^2V(^,g^ + (l-x)D^)V^2-(^,9^)- (3.15) 



Using Eq. ()3.5|) we get. 
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/ t4 — ^[(^^^ + - - ^2(1 - xi] 



+ {iD^ + D'^){1- x)xi{l- X2Y]. (3.16) 

The g-*" integration can be performed either using the Feynman parameter method or by 
using 

1 1 

d(3(3''-^e-^^. (3.17) 



A'' T{k) Jo 



Here we use the latter method and we get, 

d q — — = -ITT da d(3- — ——{l-x)Dye ' ^ 

L1L2 Jo Jo [a + p)^ 

where Dy = — iD^. We introduce the variables A and y defined as, 

\ = a + p, a = y\, P = {l-y)\. (3.19) 
The above integral can be written as, 

f^2^±^Jjn ^ f'dy{l-y){l-x)D^ r^;^^-A,(l-,)(l-.)^Dig-A,n^[,(l-x.)^+(l-,)(l-.0^] 

J L1L2 Jo Jo 

. /-i {l-y){l-x)D^ 
~ ''""Jo ^y{l -y){l- xfDl + m''[y{l - X2Y + (1 - y){l - x^f] ' ^^'^^^ 

The other integral can be done in a similar way :, 



L1L2 



1 1 
'0 "^^1/(1 -y)il- xfDl + m2[|/(l - X2Y + {l-y){l- xif] ^^'^^^ 



Substituting in Eq. ()3.18p we obtain. 



-unm 



(2vr)3 y'(l-Xi)(l-X2) 



^ J - X2f - 3:2(1 - xif]{l -y)- - Xif 

y{l-y){l-xyDl + m^[y{l-X2y + {l-y){l-xi)^y ^' ' 

In the forward limit, ^ = 0, Xi = X2 = = and we obtain, using Eq. ()2.5|1 

E,{x, = - / dy— )^ > ^ = -X, 3.23 

which gives the Schwinger value for the anomalous magnetic moment of an electron in QED 

f' E,dx = F2{0) = ^. (3.24) 
JO 2n 



IV. GAUGE BOSON GPDS 

We calculate the helicity non-fiip part of the gauge boson matrix element F^f for the same 
state as before. This gives information on the gauge boson GPDs Hg and Eg. Contribution 
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comes from the two-body wave function, which has one fermion and one gauge boson as 
constituents. This can be written as an overlap 



Tp+g 



j d\^r2+{- 



The g-*- integral is divergent. The above can be calculated using ()3.5|) : 
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A2 



— m^(l — Xi)^/ + (1 + X1X2) T^logl- 



l-xyDl + m'^{l-X2y' 

+ 2771=^(1 -X2)'(l-Xi)2/}, 



771^(1 — Xi)"^ 

(1 - x)^^! + m2[(l - X2)' - (1 - xi)2])/ 

(4.2) 



where / = / L2 = (g^ + (1 - a;)^-^)^ + m'^{l - ^2)^ and Li = {q^y + 7n^{l - Xi] 



2. 



Xi 



l-x 



1-x 



1-5' "^2 1+^- 



Like the quark case, it is again interesting to look at the forward limit of the above 
expression. We get 



FZix,0,0) = ^log 



a , A^l + (l-x) 



X 



(4.3) 



2ti ^/i2 

Here we have neglected the electron mass for simplicity. F^^(x, 0, 0) gives the (unpolarized) 
gluon distribution in the nucleon and the above expression gives the splitting function [l4 1 

1 + (1-X)2 



X 



(4.4) 



In the off-forward case, the splitting functions can be found using the same approach 

We now calculate the hehcity flip gauge boson GPD Eg given in ()2.2|1 . for the same state. 
The matrix element is given by, 

1 



SttxP^ 



dy-e^^^^y-^{P'+ I F+°(-^)F+«(^) | P-). 



(4.5) 



Contribution comes from the two particle sector. As in the fermion case, the mass cannot 
be neglected here. This can be written as, 

.1 — X 



,q^ + {l-x)D^). 



(4.6) 
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Using Eq. ()3.5|) we get, 



r,_ = - — —itm 



rfV^[xi(l-X2r-X2(l-Xi)2] 



-{\ - x)x2(\ - x^)HD^ f p^. (4.7) 

J L1L2 

The integration can be performed in a similar way as for the fermions and we get, 

F+f = ^my^l^Di^il - X) 
Ztt X 

1 , [xi(l - X2f - X2(l - xi)^](l -y)+ X2(l - x,f 
^i/(l-y)(l-a;)2Z}i + m2[y(l-a;i)2 + (l-|/Vl-x2)2]- 

In the forward limit, Xi = X2 = 1 — x and we get, in agreement with []Jj|, 

Eg{l-x,0) = -- ^ ^ ' , (4.9) 

here E'g is given by Eq. ()2.5j) . Note that x in the forward case is the momentum fraction of 
the gauge boson. The second moment of -E'q_g(a;, 0), / dxxEg^g{x,0) gives in units of ^ by 
how much the transverse center of momentum of the parton q, g is shifted away from the 
origin in the transversely polarized state. When summed over all partons, the transverse 
center of momentum would still be at the origin. Indeed it is easy to check for a dressed 
electron jl^ 

dxxEq{x, 0) + dx{l - x)Eg{x, 0) = 0, (4.10) 

which is due to the fact that the anomalous gravitomagnetic moment of the electron has to 
vanish Note that in the second term, (1 — x) is the momentum fraction of the gauge 
boson. 



V. GPDS IN THE IMPACT PARAMETER SPACE 

Fourier transform of the GPDs with respect to the transverse momentum transfer A"*- 
brings them to the impact parameter space. When the longitudinal momentum transfer ^ = 
0, this gives the density of partons with longitudinal momentum fraction x and transverse 
distance b from the center of the proton. For non-zero ^, the Fourier transforms are defined 
as ^, 

nn^^^X) = J^e-^^''^Ftr{x,^,D^) 

= ^(^^)'^o(|/^||M)(^^.,.-Y3^^.,.) (5.1) 
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-Tn—l "i^^fMlOW^nLJE,,,, (5.2, 

where Jo and Ji are Bessel functions and b-^ is called the impact parameter. In order to 
avoid infinities in the intermediate steps, we take a wave packet state 

/^0(p^)b^p^A), (5.3) 

which have a definite plus momentum. Following (3I we take a Gaussian wavepacket, 

cP{p^) = G{p^,a') (5.4) 

where 

G{p^,a^) = e^&. (5.5) 

a gives the width of the wave packet. It is the accuracy to which one can localize information 
in the impact parameter space The states centered around 60 with an accuracy a are 



normalized as, 



6'^ y I p+, 6^ A) = ^^4g(6'^ - 6^, 2a')5ip-^ - p'^). (5.6) 



Fourier transform of the matrix elements in Eq. ()5.1|) and ()5.2|) with the Gaussian wave 
packet probe partons in the nucleon at transverse position h but when the initial and final 
state protons are centered around but shifted from each other by an amount of the order 
of I fo"*" |. In our case, they probe a bare electron or a photon in an electron dressed with 
a photon. 

It is interesting to look at the qualitative behavior of the helicity-fiip GPDs Eq and Eg 
in the impact parameter space. The contribution in this case comes purely from the two- 
body sector, that is it involves the wave function il)2 of the relativistic spin-1/2 system. The 
overlap is given in terms of the light cone wave functions whose orbital angular momentum 
differ by AL^ = ±1 [3]. The scale dependence, as mentioned before, is suppressed here, 
unlike Hq and Hg. We use the notation 



£q,g{x,^M) = ^ r d{D^fM\ D\\h \)Eq^g- (5.7) 



00 
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C(^'^'^^) = i^X diD^rJii\D\\b\y-^E,,, (5.8) 

which contribute to Eq. (|5.1|) and (|5.2p respectively. Fig. 1(a) shows Sg vs b-^ for fixed 
^ = 0.1 and three different values of x > ^. We have plotted for positive b-^. The functions 
are symmetrical in | fe-*- |. For all numerical studies, we have taken a Gaussian wave packet. 
Sq is positive and has a maximum for | 6-*- |= and decreases smoothly with increasing 
I 6-*- |. For fixed \ b-^ |, is higher in magnitude for higher x. The qualitative behaviour 
for nonzero f is the same as f = 0. We have taken the normalization to be = 1 and 
m = 0.5 MeV. Fig. 1 (b) shows £q vs ^ for fixed b-^ = O.lMeV"^ and three different values 
of X > Sq is again a smooth function of ^ and increases as ^ increases for given 6-*- and x. 
For fixed ^, Sg increases as x increases. We have plotted Sq in Fig. 1 (c) as a function of b^ 
for ,^ = 0.1 and three different values of x. The rise near 6-*- = is much more sharp here 
than in Fig. 1(a). 

In Fig. 2 (a) we have shown Sg{x,^,b^) vs b-^ for fixed ^ = 0.1 and for three different 
values of X. Sg is negative for positive fe"*". It has a negative maximum at fe"*" = and 
smoothly decreases in magnitude as b-^ increases. Again, the qualitative behaviour is the 
same as = 0. For a given 6-*-, Sg decreases in magnitude for increasing x. Fig. 2 (b) 
shows Sg as a function of x for fixed b-^ = O.lMeV"^ and three different values of ^ < x. Sg 
vanishes at x = 1. For ^ = 0, it also vanishes at x = 0. For C, > 0, Sg increases in magnitude 
as ^ increases. The curves cannot be continued for x < ^ as there will be contribution 
from the higher Fock components in this region. Fig. 2 (c) shows Sg vs b^ for ^ = 0.1 
and three different values of x. The qualitative behaviour is the same as in (a), again the 
rise near 6-*- = is much sharp like the fermion case. For large b^, S^ behaves in the same 
way independent of x. For completeness, in Fig. 3 (a) and 3 (b) we show Hq{x,0,t) and 
Hg{x, 0, t) in the impact parameter space, Hqi^x, b^) and Hg^x, b^) respectively for a definite 
scale A = 5 GeV. The width of the Gaussian is a = 0.1. Both of them are smooth functions 
of 6-*", increases as | 6^ | decreases. We have omitted the very small b^ region in order to 
show the resolution of the different curves at higher b^. For a given 6^ 'Hq{x, b^) increases 
as X becomes closer to 1 and at x 1 it becomes a delta function lOj which can be seen 
analytically. T-Lg{x,b-^) on the other hand, decreases in magnitude for a given 6-*- as x goes 
closer to 1. 
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VI. INEQUALITIES 



GPDs in the impact parameter space obey certain inequalities, which impose severe con- 
straints on phenomenological models of GPDs. The most general forms of these inequalities 
were derived in from positivity constraints. The spin flip GPDs Sn r,(x, b-^) defined as 
the Fourier transform of Eq^g(x,b^) for ^ = obey two inequalities given in 19|, both of 
them can be shown to hold for a dressed electron state. For non-zero C,, a general inequality 
can be derived 0, 

(1 - er\ I'Ai^, b^) f < 0, Y^)nt(^, o, (e.i) 

for ^ < X < 1 and for any combinations of helicities A', A. For A' ^ A, the inequality is 
easy to prove, as there are large logarithmic contribution to from the scale dependent 
part. For A' = A, the inequality is non-trivial and can be verified numerically. Fig. 4 shows 
the Ihs and rhs of Eq. ()6.1|) vs for the fermions for two different choices of the scale 
A = 7 GeV and A = 5 GeV. The scale dependence comes entirely from Hq. In the plot 
xi = ^2 = fr|, bi = jz^ and 62 = We have taken ^ = 0.1 and x = 0.5 and cr = 0.1. 

VII. SUMMARY AND DISCUSSIONS 

In this work, we have investigated the GPDs for an effective spin-1/2 composite rela- 



li that 



tivistic system, namely for an electron dressed with a photon in QED. It is known 
the light-cone wave function of this two-body state gives a template for the effective quark 
spin-one diquark structure of the proton light cone wave function, which provides the phe- 
nomenological relevance of our study. The GPDs are expressed as overlaps of the light-cone 
wave functions, which in this case are known order by order in perturbation theory. We 
keep both the skewedness ^ and the transverse momentum transfer A"*" non-zero, which is 
relevant for deeply virtual Compton scattering experiments to probe the GPDs. Fourier 
transform with respect to the transverse momentum transfer brings the GPDs to the im- 
pact parameter space. We showed that the GPDs for the effective state obey the necessary 
inequahties and investigated the qualitative behaviour of the fermion and the gauge boson 
GPDs in the impact parameter space. 
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FIG. 1: (a) £g vs b-^ for C = 0.1 , (b) £g vs ^ for = OAMeV'^, (c) £^ vs b-^ for ^ = 0.1 and three 
different values of x. We have taken a = 0.001. 
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FIG. 2: (a) £g vs ft"*" for ^ = 0.1 and three different values of x. (b) £g vs x for = O.lMeV^^ 
and three different values of i^. (c) £g vs with ^ = 0.1 and three different values of x. We have 
taken a = 0.001. 
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FIG. 3: (a) 7iq{x,b-^) vs 6-*-, (b) 7ig{x,b-^) vs. for three different values of x. A is given in 
MeV. 
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FIG. 4: Inequality for /++. A is given in MeV. Here xi = bi = j^, X2 = j^, 62 = 
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